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OPTIMUM  AND  SUB-OPTIMUM  TECHNIQUES  FOR  EXTRACTION 
OF  INPUT  SIGNALS  FROM  THE  OUTPUTS  OF  NOISY 
TIME-VARIABLE  MULTIPATH  CHANNELS 


By 


F .  Bryn 

ABSTRACT 


In  underwater  communication  and  echo-ranging  the  transmitted 
signal  will  generally  reach  the  receiver  via  several  distinct 
paths.  The  resultant  signal  is  referred  to  as  a  multipath  signal. 
The  first  part  of  this  report  derives  the  integral  equation  that 
defines  the  optimum  filter  for  extracting  the  transmitted  signal 
from  a  noise-corrupted  multipath  signal.  The  solution  of  this 
equation  is  not  considered  to  be  easily  accessible.  The  second 
part  of  the  report  describes  a  practical  technique  for  deriving 
near-optimum  extraction  filters  applicable  to  echo-ranging 
situations  involving  reflections  from  a  smooth  ocean  surface. 
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INTRODUCTION 


Consider  the  randomly  time-varying  multipath  transmission  channel 
depicted  in  Fig.  1.  A  transmitted  signal  s(t)  gives  rise  to  a 
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noise-corrupted  multipath  signal 


y(t) 


N 

=  £ 
i=  1 


r 

hi(T,t)  s ( t —  t)  dT  +  n(t) 

0 


(Eq.  1) 


in  which  h  (T,t)  is  the  response  of  the  i'th  channel,  at  time  t, 
to  a  unit  impulse  applied  t  seconds  earlier. 


The  present  report  investigates  the  problem  of  defining  the  optimum 
linear  filter  for  extracting  the  transmitted  signal  s(t)  from  the 
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received  waveform  y(t) •  It  was  considered  to  be  of  particular 
interest  to  study  extraction  techniques  applicable  to  multipath 
explosive  echo-ranging  situations  in  which  it  is  desirable  to 
isolate  the  response  of  the  target  to  the  explosive  pulse. 

The  similarity  between  this  problem  and  the  pure  transmission 
problem  presented  above  is  best  understood  by  considering  Fig.  2a. 


An  explosive  type  pulse,  in  the  figure  denoted  by  6(t),  is 
transmitted  from  a  point  in  the  ocean  at  time  t  =  0.  The  pulse 
travels  via  one  direct  path  (D)  and  one  surface-reflected  path  (S) 
to  the  target  represented  by  the  two  filters  and  F^ .  The 

output  signals  from  the  filters,  being  the  responses  of  the  target 
to  the  two  incoming  pulses,  are  reflected  towards  the  receiver, 
each  via  one  direct  and  one  surface-reflected  path.  If  the  impulse 
responses  of  the  filters  are  denoted  by  Sj(t)  and  sa(t) 
respectively.  Fig.  2a  can  be  redrawn  as  shown  in  Fig.  2b.  Finally, 
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FIG.  2b  EQUIVALENT  DIAGRAM  FOR  MULTIPATH  ECHO-RANGING 

if  sl(t)  =  s2(t)  =  s(t)  the  diagram  becomes  as  shown  in  Fig.  2c. 
This  is  similar  to  Fig.  1  with  N  =  4. 


D-D 


FIG,  2c  FINAL  EQUIVALENT  DIAGRAM 

The  results  obtained  are  twofold.  Firstly,  the  integral  equation 
governing  the  form  of  the  optimum  linear  extraction  filter  for  the 
noise-corrupted,  randomly  time-varying,  multipath  channel  is 
derived;  the  solution  of  this  equation  is  not  believed  to  be  easily 
accessible.  Secondly,  a  practical,  intuitive  technique  for 
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deriving  near- optimum  extraction  filters  for  multipath,  noise- 
corrupted,  echo-ranging  situations  involving  reflections  from  a 
smooth  surface,  is  considered.  The  time-varying  element  is 
here  eliminated  and  only  the  noise  retained. 
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THE  IMPULSE  RESPONSE  OF  THE  CHANNEL 


From  Eq.  1  the  time-varying  impulse  response  of  the  channel  is 
N 

h(T,t)  =  E  hi(T,t)  ,  (Eq,  2) 

i=  1 

where  h^(i,  t)  defines  the  response  of  the  i'th  channel,  at 
time  t,  to  a  unit  impulse  applied  t  seconds  earlier.  In 
physical  systems  we  always  have 


h  ( T,  t)  =0  for  t  <  0 


The  Fourier  transform  of  h(T,t)  with  respect  to  t  is  the 
time-varying  frequency  response  H(f,  t)  of  the  channel.  Thus 


H(f ,  t) 
and 


r 

J  h( T ,  t ) 

0 


-  j  2rrf  t 
e 


dr 


h(  T,  t) 


j  H(f ,  t)  ej2TTfT  df  . 


(Eq.  3) 


In  what  follows  we  shall  assume  that  h(T,t)  and  H(f,  t)  are 
random  functions  of  time  t  and  that  they  are  members  of  stationary 
processes  h(<r,t)  and  H(f,t)  respectively. 
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THE  OPTIMUM  EXTRACTION  FILTER 


Referring  to  Eq.  1,  the  problem  considered  is  as  follows: 
Given  a  mixture 


y(t)  =  m(t)  +  n(t) 


(Eq.  4) 


of  a  multipath  signal 


« 


s ( t  -  t)  dT 


(Eq.  5) 


and  stationary  noise  n(t) ,  what  is  the  optimum,  physically-realizable, 
linear  filter  for  estimating  s(t  -  tc)  from  y(t)?  The  delay  t0 
may  be  arbitrarily  selected  by  the  filter  designer. 

It  is  assumed  that  the  transmitted  signal  s(t)  has  finite 
length  T  and  that  we  only  are  interested  in  the  filter  estimator 
output  in  the  interval  from  t0  to  (t0  +  T) . 

Let  g( t)  be  the  impulse  response  of  the  filter  estimator  and 
Z(t)  its  output,  as  shown  in  Fig.  3»  Mathematically  our  problem 


y(t)=m(t )+n(t) 


g(  t) 


Actual  output  =  z (t) 
Desired  output  =  s(t  -t0) 


FIG.  3 
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is  identical  to  that  of  determining  g( t)  such  that  the 
mean  square  error: 


g 


(Eq.  6) 


is  a  minimum.  Note  that 


z(t)  = 


r 

Jo 


g(  t)  •  y(t  -  t)  dx  . 


The  minimization  may  be  done  by  variational  techniques,  i.e. 
one  introduces  a  new  filter  estimator  with  impulse  response 

k  (  t)  =  g(  t)  +  e  •  t)(  t) 

and  requires  that  the  derivative  of  the  corresponding  mean  square 
error  with  respect  to  e  be  zero  for  e  =  0  for  all  functions  1"|(t). 

It  is  assumed  that  the  properties  of  T|(t)  are  those  of  a 

physically  realizable  filter  response. 

The  minimization  process  is  shown  formally  in  Appendix  A.  The 
result  is  the  integral  equation: 


dT  -  CO  (t„  , v)  =  0  v  >  0 
^m . s  °  - 


(Eq.  7) 
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where 


t0+T 


The  function  m(t)  is  given  by  Eq.  5*  Figures  4a,  b,  c,  and  d 
present  some  simplified  illustrations  of  the  signals  involved  in 
Eq.  7. 

If  t0  is  chosen  sufficiently  large,  all  terms  on  the  lefthand 
side  of  Eq.  7  will  approach  zero  for  v<0.  It  is  then  possible 
to  apply  a  Fourier  transform  with  respect  to  the  variable  v. 
Introducing 


G(f  ) 


NU) 


J 


dT  •  dv  , 


and 
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s(t) 


m(t) 


/ 


\ 


XA7 


Output  interval 
/ - A  \ 


H - \ — ►  t 

t0  t0+T 


FIG.  4a  TRANSMITTED  SIGNAL  s(t)  AND  MULTIPATH  OUTPUT  m(t) 
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FIG.  4b  THE  CORRELATION  FUNCTION  cp  (t,  v) 

mm 
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T 


FIG.  4c  THE  CORRELATION  FUNCTION  Cp  (t-  v) 

nn 
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m.s 
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o  > 


V) 


Jy  Jy -A- 


V 


FIG.  4d  THE  CORRELATION  FUNCTION cp  (t a,v) 

m.  s  ° 
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Eq.  7  may  be  transformed  to  yield 


p+» 


G(f )  •  $  (f ,£)  df  4-  T  •  G( £)  • 

m  o  m 


N(^) 


$ 

m.s 


(t 


o  > 


u  . 


Introducing  finally 


F 

m  m 


(f,0 


$  m  (f ,5)/T-N(5) 
mm 


and 


F 

m.s 


(t0  f 


one  obtains 

p+» 

G(f)*  F  (f,5)  df  +  G(5)  =  F  (tc,0 
.  I  m  .  m  in . s 

w«oo 


(Eq.  8) 


This  will  be  recognized  as  an  inhomogeneous  integral  equation 
of  the  second  kind. 
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A  PRACTICAL  EXTRACTION  TECHNIQUE 


The  theoretical  treatment  of  the  extraction  problem  presented 

in  Ch.  2  can  hardly  be  considered  very  practical  for  echo-ranging 

situations,  because  the  signal  and  path  information  required  for 

the  functions  F  (f,  £,)  and  F  (t0,  £)  of  Eq .  8  will  rarely  be 

m  .  m  m  o  s 

available.  In  this  chapter  we  shall  attempt  to  solve  the  extraction 
problem  for  an  echo-ranging  situation  involving  time-invariant 
transmission  paths  by  means  of  a  more  practical  intuitive  approach. 

The  situation  to  be  considered  is  similar  to  the  one  depicted  in  Fig. 
2c o  It  is  assumed  that  the  paths  are  time-invariant,  and  that, 
except  for  the  path  delays  which  may  be  different,  they  have 
identical  and  ideal  transmission  characteristics.  Such  conditions 
may  be  encountered  when  the  ocean  surface  is  very  calm.  Furthermore, 
it  is  assumed  that  the  points  of  transmission  and  reception  are 
the  same.  Paths  D-S  and  S-D  will  thus  coincide  and  the  impulse 
response  of  the  path  combination  between  points  a  and  b  in  Fig.  2c 
is  as  shown  in  Fig.  5a .  The  absolute  value  of  the  Fourier 
transform  H(f)  of  this  impulse  response  is  as  shown  in  Fig.  5b. 

It  is  easily  shown  that 


-j  2nfta 

H(f)  =  -2(1  -cos  2tt  •  f  •  At )  •  e 


( Eq.  9a) 


where  At  is  the  separation  between  successive  paths  and  t3  is 
the  delay  associated  with  the  D-S/S-D  paths.  Note  that  the 
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D-D 


D-S 

S-D 


S-S 


FIG.  5a  IMPULSE  RESPONSE  OF  MULTIPATH  CHANNEL 


FIG.  5b  FREQUENCY  TRANSFER  FUNCTION  OF  MULTIPATH  CHANNEL 
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pattern  of  Fig.  5a  is  symmetric  about  the  D-S/S-D  paths  and  that 
the  responses  due  to  these  are  inverted  due  to  one  surface 

-j  ^  ^ 

reflection.  The  time  delay  expressed  by  the  term  -e  2  in 

H(f)  will  be  ignored,  since  this  is  of  no  importance  for  what 
follows.  Thus  we  shall  write 


H(f) 


2(1-  cos  2tt  f  •  At ) 


(Eq.  9b) 


The  idea  behind  the  extraction  technique  will  be  described  by 

reference  to  Fig.  6a.  A  signal  s(t)  equal  to  the  impulse 

response  of  the  target  is  the  input  to  a  linear  filter  with  impulse 

response  equal  to  the  multipath  structure  shown  in  Fig.  5a.  Gaussian 

noise  n(t)  is  added  to  the  multipath  signal  and  the  mixture 

fed  to  the  inverse  filter,  the  structure  of  which  is  to  be 

determined.  The  output  of  this  filter  can  be  considered  to  be  the 

sum  of  three  components.  The  first  is  a  desired  signal  k°s(t). 

The  second  is  undesired  signal  components  s  (t).  The  third 

m 

is  the  output  nQ(t)  due  to  the  noise  n(t).  The  three 
components  are  indicated  in  Fig.  6b.  The  design  criterion  for  the 
inverse  filter  is  that  it  should  maximize  the  ratio  between  the 
energy  of  the  desired  signal  and  the  energy,  in  the  same  interval, 
of  the  undesired  signal  components  and  noise.  Thus  one  wants  to 


maximize  the  ratio 


t0+T 


k/ 


r 

4. 


sa(t)  dt 


t  +T 


E(j  |>m<t)+n0(t)]  dt) 
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Multipath 

Structure 


s(t) 

J 

*\ 

-<+> 


Inverse 

Filter 


Desired  signal 
k  s(t ) 

Undesired  signal 

components  s  (t) 
m 


_*.Noise  n  (t) 


n(t) 


FIG.  Sa  EQUIVALENT  DIAGRAM  OF  MULTIPATH/ INVERSE-FILTER  COMBINATION 


k*s(t) 


t  +T 

O 


s  (t) 


■  m  i 


n«(t) 


-n-J-wWV 


FIG.  6b  THE  COMPONENTS  AT  THE  OUTPUT  OF  THE  INVERSE  FILTER 


where  the  statistical  averaging  indicated  by  E  ^  J  applies  to 

the  noise  n0(t).  Since  the  noise  will  be  statistically 

independent  of  s  (t)  and  will  be  assumed  to  have  zero  mean 

m 

value,  one  can  write 

k3  .  P 

r  “  777  2  -  (E<>-  101 

T  I  dt  + 

where 

1 0+T 

Ps  =  T  |  sS(t)  dt  • 

We  shall  first  consider  the  frequency  function  H.(t)  of  the 
ideal  inverse  filter  when  n(t)  =  0,  that  is  for  the  noiseless 
situation.  The  form  of  this  filter  is 

H.(f)  =  1/H(t)  , 

with  H(f)  given  by  Eq.  9  as 

H(f)  =  2(  1  -  cos  2tt  f  *At)  . 

Thus 

iMf)  =  1/2(1-  cos  2rr  f  •  At )  . 
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The  functions  H(f)  and  H^f)  are  shown  in  Fig.  7.  The  product 
of  the  two  is  of  course  equal  to  unity  for  all  f  and  is  the  desired 
total  response  yielding  an  output  equal  to  s(t) . 


FIG.  7  THE  IDEAL  INVERSE  FILTER  n(t)  -  0 

Consider  next  the  more  practical  situation  with  n(t)  0.  We 
shall  assume  that  n(t)  has  uniform  power  spectral  density. 
Referring  to  Fig.  7  it  will  be  apparent  that  the  inverse  filter, 
with  its  large  amplitudes  at  frequencies  that  are  multiples  of 
l/At,  will  introduce  noise  of  very  high  intensity  at  the  output. 
It  is  therefore  necessary  to  modify  the  inverse  filter  in  a 
direction  that  will  tend  to  maximize  the  ratio  r  of  Eq.  10. 
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Figure  8  illustrates  one  possible  modification.  It  is  seen  that 
H^(f)  is  periodic  with  period  l/At . 


FIG.  e  MODIFIED  INVERSE  FILTER  FUNCTION 


For  f<  l/2At  we  have 

(l/2  ( 1  -  cos  2n  f . At)  for  “  c  |  f | <  l/2At 

(  2 

Hi(f)  =  )  (Eq.  11) 

(  0  for  |f  |  <  “ 
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The  problem  is  now  to  determine  b  such  that  the  ratio  r  of 
Eq .  10  is  maximum.  It  should  be  noted  that  this  is  only  one 

possible  form  of  the  inverse  filter.  Other  modifications  involving 
functions  with  less  abrupt  changes  may  be  more  useful. 

Figure  8  also  shows  the  total  frequency  response  H(f)#H^(f)  This 
may  be  decomposed  into  the  two  components  shown  in  Fig.  9. 


Thus 


FIG.  9  DECOMPOSITION  OF  H(f)*H  (f) 


H(f  )  •  H  (f )  =  1  -  H  (f ) 

i  m 


(Eq.  12) 


For  practical  purposes  it  is  necessary  to  impose  bandwidth 
limitations  on  the  inverse  filter  or  the  receiving  circuits 
preceding  it.  We  shall  assume  that  the  signal  s(t)  contains  no 
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frequency  components  at  frequencies  greater  than  f  =  B,  where 
B  >>  1/At  , 

and  that  the  receiving  circuit  preceding  the  inverse  filter  has 
the  characteristics  of  an  ideal  lowpass  filter  of  bandwidth  B. 

If  the  input  noise  power  in  this  frequency  band  is  it  is  easily 

shown  with  the  aid  of  Appendix  B  that  the  signal-to-noise  ratio  r 
of  Eq.  10  becomes 


oe  m 


(Eq.  13) 


72  72  C  .  .  •  a  .  •  a  .  +  —  |  tan  ^  (  1-b  •  At )  +  —  tan  ^  (  1-b  •  At ) 

,  ,  .  .  ij  i  J  8rr  L  2  3  2  J 

i=l  j=l 


where 


t  o+T 


Ps  “  T  {  sS(t)  dt 
t 


t.+T 


C_=  ^  f  |s(t-t0-i*  At)+s(t-tc+i.  At)  l|s(t-t0-j  *  At)  +s(t-t0+j  .At)  J 
t 


and 


a0>  a£  are  zero*th  and  i*th  coefficient  respectively  in  the 

Fourier  Series  expansion  of  H  (f).  (See  Appendix  B.) 

m 

They  depend  on  b. 
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The  signal-to-noise  ratio  as  given  by  Eq .  13  must  be  maximized 

with  respect  to  b.  This  process  can  generally  be  greatly 
simplified  when  one  further  assumption  can  be  made.  Thus  if 


C.  . 
ij 


(  1 

( 

( 

(  0 

( 


for  i  =  j 

for  i  ^  j 


(Eq.  14) 


we  can  write 


m  oe 

z  z 

i=l  j=l 


c  .  . 

IJ 


a .  •  a 
1  3 


According  to  Appendix  B 

ao  =  b*At 
and 

00 

Z  a. 3  =  b . At (  1  -  b • At )  . 

i=  1  1 


Therefore 


—  (  1  -  b .  At ) 
N 


r  = 


'  1 

~  *b • At ( 1-b . At )+  [tan  f(  l-b*At)+^tan3  ^  (  1  _b«At)  J 

N 


(Eq.  15) 
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Equation  15  can  be  differentiated  with  respect  to  b  in  order  to 

determine  the  value  b  =  b0  yielding  the  maximum  value  of  r. 

The  value  of  b.  will  clearly  depend  on  the  ratio  P  /P.T,  as 

s  N 

will  also  r 

max 

In  order  that  the  assumption  contained  in  Eq.  14  shall  be  valid 

it  is  necessary  for  the  signal  bandwidth  B  to  be  much  greater 

s 

than  l/At  and  for  the  signal  length  T  to  be  much  greater 

than  l/b0 .  The  condition  T  >>  l/bc  ensures  that  the  major  part 

of  the  components  of  s  (t)  of  Figs.  5  and  6  will  fall  within  the 

m 

time-interval  ^t0,(t0  +  T)  j.  will  then  involve  auto-correlation 

functions  of  s(t)  with  arguments  (i  ±j)  At.  Since  the  minimum 

value  of  (i  ±j)*At,  i  5/ j  ,  is  At,  the  condition  B  >>  l/At 

s 

implies  that  these  auto-correlation  functions  are  zero  except  when 
i  =  3  • 
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4.  COMMENTS 


An  extraction  technique  following  the  very  general  theoretical 
prescriptions  of  Ch.  2  is  not  considered  very  practical  due  to  the 
amount  of  a  priori  knowledge  required  about  channel  and  signal 
(target)  statistics  and  also  due  to  difficulties  in  solving  the 
integral  equation  involved. 

The  less  general  technique  outlined  in  Ch.  3  appears  to  be 

of  greater  practical  value.  However,  further  studies  are  required 

before  the  true  value  of  this  approach  can  be  properly  assessed. 

These  studies  should  relate  to  the  accuracies  with  which  the 

quantities  P  /P  and  At  of  Eqs .  13  or  15  can  be  estimated  and 

S  N 

to  the  influence  of  estimation  errors  on  the  mean  square  error 
of  the  extracted  signal. 
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APPENDIX  A 


THE  MINIMIZATION  PROCESS 


It  is  desired  to  find  the  physically-realizable  filter  impulse 
response  g(iO  that  minimizes  the  expression 
tc+T 

=  E  (j1  £  J  g(T)*y(t-T)  dT-s(t-te)j  dt  ^  Eq.  A.l 


Inserting 


k(  t)  =  s(  T )  +  e*n(T) 


Eq  .  A  .  2 


for  g( t)  one  obtains 


t  o+T 


-  E( 


r  r  r 
j  j  j 


k( t) *k( v) *y( t-T) *y( t-v)  dv*dT*dt 


te+T 


2  £  s(t-tj 


t  _+T 

O 


f  k(T).y(t-T)  dT*dt  +  F  sa(t-t  )  dt^ 


Eq .A  •  3 


25 


After  substitution  for  h(x)  from  Eq.(A.2)  it  is  desired  to 
evaluate  the  differential  S  p  /de  at  e  =  0. 

rC 


2 


]>>[«>•  >+T 


E(y(t-x)  *y(t-v)) 


dt  *dx 


tc+T 

j  s(t-t0)  'E^y(t-v)^  dt  dv 


It  is  required  that  this  be  zero  for  all  T](v) .  This  is  obtained 
when  the  quantity  inside  the  square  parenthesis  is  zero.  Thus 


t0+T 


Jo 


g( t)  *  j£  E  ^  * y ( t- t)  *y ( t-v)^  dt  dT 


t  +T 

O 

I  s(t-t")‘E(y(t-v)) 


dt 


=  0 


for  v  >  0 


Eq .  A. 4 


Since  y(t)  is  the  sum  of  a  multipath  signal  m(t)  and 
stationary  noise  n(t)  ,  we  can  introduce  the  functions 

1 0+T 

CP  ( T 
m .  m 


,v)  = 


j  E  ^ m(t-x) *m( t-v)j  dt 
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cpn  n(r-v)  =  E  ^  n(t-r)  *n(t-v)^ 
t0+T 

CPm  s(te,v)  =  j'  s(t-t0)°  E  ^m(t-v)^  dt 

to 


into  Eq.(A.4).  The  result  is 


r 

Jo 


g(T)  r cpm  (tjv) 

L  m.m 


+  T*cp„  (t-v)"1 
n  o  n 


dT 


Cp 

m*s 


(t 


o  y 


v) 


0 


for  v  >  0 
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APPENDIX  B 


THE  OUTPUT  SIGNAL-TO-NOISE  RATIO 


The  total  response  of  the  channel-filter  combination  was  given  by 
Eq *  12  of  Ch .  3  as 


H(f ) *H . (f )  =  1  -  H  (f ) 

l  m 


Eq  .  B  .  1 


We  shall  study  this  response  function  in  the  time  domain  and 

separate  it  into  two  components,  one  giving  rise  to  the  desired 

signal  and  the  other  to  undesired  signal  components  s  (t) • 

m 


Consider  first  the  function  H  (f)  depicted  in  Fig.  B.l,  It  is 

m 

periodic,  with  period  l/At,  and  can  be  expressed  by  the  Fourier 
series 


H  (f) 
m 


a0  +  2  £ 
i=  1 


a  •  cos  2TT.i*At  f 

l 


It  can  be  readily 


shown  that 


1/2  At 

P 

a  =  At  H  (f)  df 

J  m 

-1/2- At 


b  •  At 


Eq  .  B.  2 


Eq  .  B . 3 


and  that 

l/2At 

n  a  * 

At  |H  (f)  I  df  =  a*  +  2  £  a?  . 

J  m  n=1  1 

-  1/2 • At 
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The  latter  yields 


b  •  At  =  a  3  +  2  72  a  .3 

i=  1  1 


or 


2  E  a.3  =  b’At^l  -b*At)  .  Eq,  B.4 

i=  1  1 

We  return  now  to  Eq.  B.l.  Inserting  for  H^(f)  from  Eq.  B.2 
gives 

as 

H(f)*H.(f)  =  (l-a0)  -2  E  a  cos  2tt  •  i •  At •  f 

i=  1 


Performing  a  Fourier  transform,  we  obtain  the  impulse  response 
h^(t)  of  the  channel-filter  combination.  This  function  is 
illustrated  in  Fig.  B.l  and  is 


h  (t)  =  (l-a#)6(t) 

Lj 


-  Z/  a_^  ^  6(  t-i  •  At  )  +  6( t+i  •  At )  J 


E(  B.5 


Here  the  term  (1  -  a0 )  6(t)  carries  the  desired  signal  component. 
The  remaining  terms  carry  undesired  signals  s  (t) . 


s 
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In  the  derivation  of  Eq.  B.5  ,  time-delays  in  the  channel  and 
filter  response  functions  have  purposely  been  omitted  on  the 
grounds  that  they  do  not  affect  the  magnitudes  of  the  quantities 
considered  in  Eq,  10.  A  real  situation  involves  a  channel  delay, 
which  in  Eq  .  9a  is  denoted  by  t2 ,  and  a  filter  delay  t  which 
is  necessary  if  the  inverse  filter  is  to  be  physically  realizable. 
The  delay  t^  must  be  approximately  half  the  extension  of  the 
Fourier  transform  of  H^(f)  in  the  time  domain.  Writing 

Eq ,  B.5  should  be  written 

00 

h  (t)  =  (  l-a0) 6(t-t0)~  E  ai[6(t-t0-i*At)  +  6(t-t0+i*At)]  Eq.  B.6 

E  i=  1 

If  we  apply  a  signal  s(t)  to  the  channel  input  at  time  t  «*  0, 
the  filter  output  signal  becomes 

u  ( t )  «=  (  1  -  a  )  s  ( t — 1 0  )  +  s  (t)  , 

v  v  m 

where 

00 

s^(t)  =  -  £  a^j^s(t-t0-i*  At)  +  s(t-tc+i .  At )  J  . 

i—  1 
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We  are  now  able  to  extract  the  quantities  required  for  Eq.  10. 
Thus 


k3  =  (  1  -  a0  )  3 


Eq .  B . 


1 

T 


1 0+T 


J 


s3(t)dt  =  2 
m  i=  1 


eo 

2 


c  .  . *a  *a  .  , 

iJ  i  3 


where 


te+T 


C_  =  ^  j*  j^s(t-t0-i*  At )+s(t-t0+i  •  At )  J|s ( t-t e-j  •  At )+s( t-t e+j  .  At )  J 
t. 

(B.8) 


The  final  information  needed  for  Eq*  10  is  the  evaluation  of  the 
noise  term  E  (  n0a(t)  J  *  We  can  write  this  as 


nD  2 

|H.(f)|  df  . 

0 


Since  H^(f)  is  periodic,  with  period  l/At  and  B>>l/At,  we  can 
write 

l/2 • At 

p  3 

E  (ne3(t))  =  2-  -f  •  B-At  j  |H.(f)  |  df  . 


dt 
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1/2  (  1  -  cos  2tt  *f  «At) 


Inserting  H^(f)  = 


and  performing  the  integration,  yields 


E 


8  -  tt 


l^tan -(  1  _  b.  At)  +  “tan3  “(l-b’At)J  . 


Eq .  B .  9 
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